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Abstract 


Hadron masses are shown to be separable in QCD into contributions of 
quark and gluon kinetic and potential energies, quark masses, and the trace 
anomaly. The separation is based on a study of the structure of the QCD 
energy-momentum tensor and its matrix elements in hadron states. The pa¬ 
per contains two parts. In the first part, a detailed discussion of the renor¬ 
malization properties of the energy-momentum tensor is given. In the second 
part, a mass separation formula is derived and then applied to the nucleon, 
pion, and the QCD vacuum. Implications of the results on hadron structure 
and non-perturbative QCD dynamics are discussed. 
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I. INTRODUCTION 


For any physical system, a good knowledge about its mass structure is helpful in un¬ 
derstanding the underlying dynamics of the system. The ground state energy of nuclei as 
a function of the mass number reveals the nuclear shell structure and the saturation prop¬ 
erty of nuclear force [|l|. An organized pattern of hadron masses led to the speculation of 
the quark substructure of hadrons [|| . Likewise, any knowledge about structure of hadron 
masses in terms of its underlying constituents, quarks and gluons, can be useful to unlocking 
the physics of Quantum Chromodynamics (QCD) in strong coupling region. 

At classical level, the lagrangian of Chromodynamics is invariant under scale transfor¬ 
mation if quark masses are neglected, and thus hadron masses necessarily vanish. Radiative 
quantum effects break the scale symmetry J3| and introduce a dimensional parameter Aqcd 
through dimensional transmutation |ij. [For a classical discussion on the scale anomaly, see 
Ref. [3|]. Although Aqcd is well-determined through scale-breaking effects in high-energy 
processes 0, the physical mechanism for generating the scale at low energy is not quite 
clear. [It is encouraging, though, that the scale determined from hadron spectrum using 
lattice QCD is consistent with that determined at high energy [B[].] 

In the past, scale generation at low energy is largely understood in two seemingly-unrelated 
pictures, and so are masses of hadrons. The first picture emphasizes the aspect of color con¬ 
finement, through which quarks in hadrons are confined to a cavity of radius ~ 1 fm. A 
simple, representative model is the MIT bag, in which the mass of the nucleon is generated 
from the quark kinetic energy and vacuum pressure 0 . A dimensional analysis shows that 
the quark kinetic energy accounts for 3/4 of the mass and the vacuum pressure accounts for 
the rest 1/4. In some loose sense, the vacuum pressure models effects of long-wavelength glu¬ 
ons. The second picture for scale generation emphasizes chiral symmetry breaking. Through 
such phase transition in the QCD vacuum, light quarks acquire a constituent mass of order 
300 MeV. Masses of hadrons are then approximately the sum of the constituent masses and 
quark kinetic and potential energies || . Both pictures work quite well for hadron spectra and 
other physical observables. It is not known, however, which picture or which combination 
of the two is closer to reality. 

In a Letter article published recently [0, this author showed that an insight about the 
mass structure of hadrons can be obtained through a study of the energy-momentum tensor 
of QCD. The result is a separation of hadron masses into contributions from quark and gluon 
kinetic and potential energies, the current quark masses, and the trace anomaly. The last 
part is a direct result of scale symmetry breaking, and is analogous to the vacuum pressure 
empirically introduced in the MIT bag model. Though the scale of other contributions is 
determined by the anomaly, relative magnitudes reflect important aspects of the low-energy 
quark-gluon dynamics. The separation here is analogous to the virial theorem for a simple 
harmonic oscillator and the hydrogen atom. 

The goal of this paper is to provide more field theoretical basis for the mass separation 
and to extend its application to other hadrons besides the nucleon. The field theoretical 
discussion mainly answers questions like: Why the QCD hamiltonian is finite? Why it 
can be separated into gauge-invariant and finite pieces? How the renormalization scale 
dependence affects the separation? Why there is an extra term in the hamiltonian from the 
trace anomaly? etc. An application of the mass separation to a specific hadron requires 
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knowledge of two matrix elements: the momentum fraction of the hadron carried by quarks 
in the infinite momentum frame and the quark scalar charges. The matrix elements are 
known to a good accuracy in the nucleon and pion. However, there is no firm estimate of 
these for other hadrons. Nonetheless, a plausible assumption leads to a crude picture for 
general mass partition. 

The paper consists of two main sections. In Section II, we start with the separation of 
the QCD energy-momentum tensor into the traceless and trace parts. Then we consider 
renormalization of both parts in dimensional regularization and the covariant gauge fixing. 
We present a derivation the trace anomaly from the point of view of operator mixing. We also 
discuss structure of the mixing matrix for the traceless part of the energy-momentum tensor. 
As a by-product, we derive the dilatation current in QCD and study its anomalous Ward 
identities (Callan-Symanzik equation). For a phenomenology-oriented reader, this section 
can be skipped. In Section III, we study the hadron matrix element of the tensor and 
derive a separation of the QCD hamiltonian and hence hadron masses. For the application 
to the nucleon, we focus on discussing physical significance of the different contributions. 
The result for the pion strongly supports the concept that the pion is a collective vacuum 
excitation. Finally, comments are made about general features of the mass separation for 
other hadrons and the QCD vacuum. 


II. RENORMALIZATION OF THE QCD ENERGY-MOMENTUM TENSOR 


In this section, we discuss renormalization properties of the QCD energy-momentum 
tensor. A similar but less extensive discussion was first made by Freedman, Muzinich, and 
Weinberg |l(J for non-Abelian gauge theory. Since then, many papers relevant to the subject 
have appeared in the literature. The present discussion is not a review of the subject. In 
particular, there are still open questions that are under current debate. Rather, we will focus 
on the following three aspects that are most useful for our discussion of hadron masses. First, 
the energy-momentum tensor can be uniquely separated into the trace and traceless parts 
which are to be separately renormalized. Second, the trace part of the tensor is shown 
to be finite according to that the renormalized Green’s functions are finite functions of 
renormalized masses and couplings. The result is equivalent to a new derivation of the trace 
anomaly originally given in Refs. ||TT|,|T2f. Third, the traceless part of the tensor is shown 


to be finite using Ward identities related to space-time translational symmetry. A study of 
scale symmetry and related Ward identities (Callan-Symanzik equation) is included in the 
end of this section as a by-product of the previous discussion. 

The fundamental QCD lagrangian reads, 


£qcd = — m )-0 - ^F iw F iju , 


( 1 ) 


where if} is quark fields, carrying implicit flavor, color and Dirac indices, D M + igA M 

is a covariant derivative with A M = kF a f a being the gauge potential (Tr t a t b = 8 ab / 2), m is 
the quark mass matrix in flavor space. The gluon held strength has the usual non-Abelian 
expression, 


jpfii/a _ _ Q v A^ a gf abc A^ b A UC 


( 2 ) 
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where f abc is the structure constant of the gauge group. To simplify notations, we shall 
omit color indices if no confusion arises. All quantities without further specification are bare 
ones. 

To discuss renormalization, we work in dimensional regularization and (modified) min¬ 
imal subtraction scheme. An important virtue of this scheme is that it does not modify 
the basic QCD lagrangian except all Lorentz indices and space-time integrations are taken 
to be d dimensional. [The Pauli-Villars regularization, some momentum cut-off schemes, 
or lattice cut-off all modify the lagrangian in a more significant way.] We use perturbation 
theory to explore large momentum behavior, and thus a gauge fixing is needed. We choose 
the usual prescription for use of the covariant gauge, which introduces two extra terms in 
the lagrangian: the gauge fixing term, 

4f=-^ ( svg 2 , ( 3 ) 

and the Faddeev-Popov ghost term, 


£ gs = d^uD^uj . 

Therefore the total effective lagrangian, 

£eff £qCD T £gf T £gs i 


(4) 

(5) 


enters the rest discussion. 

The energy-momentum tensor can be derived from the fact that the action S = / d d x£ e g 


is invariant under space-time translation. The result is well-known [FLO 


T^ u = 


Aff - F ,ia F V a 

- g^C 1 d a {A a d ■ A) + 2C 1 A^d u \d ■ A) 
SS 




2 - XX A„ - Ai[y i7lv , - £-[Y,Y 


5A : 


( 6 ) 


where (/iz/) means symmetrization of the indices and D = — d + igA^. There are two 
standard methods to derive the above expression. First, one can write down the canonical 
Noether current associated with space-time translational symmetry and then improves it 

nru 


with the Belinfante procedure 


Second, one can follow Ref. 


deriving a current 

coupling to gravity. Both methods yield the same tensor for QCD, though not for the 0 4 
theory. 

The last three terms in ([|) are usually ignored in other references because they vanish 
when the equations of motion for gluon and quark fields are used. When they are inserted 
into a Green’s function, their role is to change polarizations of external lines and thus they 
are finite operators. When they are included in the tensor, the Ward identities related to 
space-time translation symmetry take a simpler form 
momentum tensor is symmetric, 


12], Without these terms, the energy 


rpfU/ _ 


(7) 


4 


















In the following discussion, we concentrate on the symmetric part of the tensor only. 

The symmetric part of the energy-momentum tensor can be split into a sum of traceless 
and trace parts. Although a splitting like this is not unique in general, we define the splitting 
in the following unambiguous way, 




( 8 ) 


where the traceless (T 7 *") and trace parts (T^) are, 


rp/J,U _ rjgflV _ _ T 

a 


OL 

a i 


GT 


(9) 


Notice that we have defined the trace consistent with dimensional regularization, under 
which tensors of different ranks do not mix. Therefore the necessary and sufficient condi¬ 
tion for the energy-momentum tensor to be finite is that the traceless and trace parts are 
separately finite. This is indeed the case as we shall see below. 


A. The Trace Part of the Energy-Momentum Tensor 

Let us first consider the trace-part of the energy momentum tensor, 


T a = ipmip + e 


Op + O 


gf 


— (2 — e)d a £ 1 A a d ■ A + a iD a u 


. ,-SS , ,_SS f e\ A SS 
~ ( 3 - - (2 - e)u— - (1 - -) A^ - 


where we have defined, 




5co 


5 A,. 


( 10 ) 


Of = -~F a0 F a0 , 

Ogi = —( 2£) -1 (<9 • A) 2 , 

1 SS 

2 ' 


Ogi — O g f Afj. 


( 11 ) 


In Refs. it was shown that the trace can be expressed in terms of a linear combination 

of a set of renormalized operators with finite coefficients. Here we would like to derive the 
same result from the operator mixing point of view [14], by extending the discussion in 
Ref. EH. 

First consider the equations-of-motion-related operators, 


O a — A, 


5S 

SA t 


= -0^ = I) —m)ij) , 
A Q 

O u = uj— = uid a D a uj . 
ou> 


( 12 ) 
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A zero momentum insertion of these operators into a renormalized Green’s function pro¬ 
duces the number of external lines. Therefore all these operators are finite and need no 
renormalization. 

The total derivative operator in the third term in Eq. © is a BRST variation of 
d a (A a u), i.e., 


<9“<5brst (A a u) = d a £ A a d ■ A + uD a u> 


(13) 


As such it has vanishing physical matrix elements at non-zero momentum HTTj| . |T7| . At zero 
momentum, the operator vanishes identically as it contains no poles. Since other dimension-4 
operators do not mix with it, it is a finite operator by itself. 

Renormalization of the rest three operators can be worked out by studying the renor¬ 
malized Green’s function generating functional, 


v,x,x) 

f D ( A R , -0 R , 7p R , uj r , u; R ^j exp i 


S + f (J^A R + r]ifj R + + x^ R + a ,R X ) 


/ D ( A R , ip R , ip R , cu R , uj R ^j exp iS 


(14) 


where the renormalized helds are related to bare fields by the multiplicative renormalization 
constants, 


-0 = Z 2 0 


R 


uj = Zu 


R 


A„ = Z,A ? 


(15) 


Since Z( J, rj, fj, x, X) is finite, its derivatives with respect to the renormalized quark masses 
m R , the gauge coupling constant g R , and the gauge fixing parameter £ R are also finite. In di¬ 
mensional regularization and minimal subtraction scheme, all the renormalization constants 
are independent of quark masses and the renormalization factors for the gauge coupling Z g 
and the quark masses Z m are independent of the gauge-fixing parameter |I 8 | . According to 
these, we find the following quantities are finite, 


O m = , 




R 


dhlZ:> 




0?= 1+9 


R 


d In Z, 


£ 


dg R j 
d In Z 2 


R 


9 '0 F +[g R 
1 


di R 

dlnZ n 


d£ R 


2 9 


R 


dg R 

d In Z 2 x 




R 


dg R 


Oib + £, 


di R 

d hi Z 3 

d In Z 

- a 


+ 2 9 


1 R dln Z 3 \Q' 


R 


2 9 


R 


dg R J 
d In Z ' 


dg R 


gf 


a 


l^ R dln Z m 

-\Jr 


(16) 


2 dg R 

From the above equations, we can solve Op and O g f in terms of the renormalized ones, 
e (Op + O g f^ = —2/3/ g R ^j O r + (—2/3/ g R + 73 ) 0 R { — 7 CG — 72 + 7 m O m , (17) 

where we have defined the anomalous dimensions, 
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72,3, m 

7 

/3 

Since all these quantities have no 1/e poles, the right hand side of Eq. ([17]) is finite. 

Inserting Eq. (|I7j) into Eq. O. we find the trace part of the energy-momentum tensor 
expressed in terms of finite, renormalized operators, 

T“ = (- 2 P/g R ) (Op + D§) + (1 + 7m ) O m - 2d a 5 BRST (A a u) 

+ 73 Og — ^1 + Oa — (2 + 7 ) — (3 + 72 ) Oq . (19) 

Two comments are in order at this point. First, composite operators defined through path- 
integral formalism have implied subtraction of their vacuum expectations. Thus they are 
said “normal ordered,” but not in the usual sense of relative to perturbative vacuum. Second, 
although individual term in the above equation depends on the renormalization scale (/x 2 ), 
but the sum does not. 

According to the Joglekar and Lee theorems UM3- the matrix elements of Oa, O^, and 
O w , and BRST exact operators vanish in a physical state. Thus as far as hadron matrix 
elements are concerned, we effectively have, 

n = (-2 (5/g R ) Of + (1 + 7m) O m • (20) 

This suggests we write the trace part of the energy-momentum tensor as, 

f = fr (7) + fc (7), ( 21 ) 

where T£“ = - (2p/g R 'j 0§ (s'"/ 4 ) and = (1 + 7 m ) O m (g^/i). 


= d 


d In Z- 


2,3 ,m 


ci/i 
d In Z 


= fi- 


dji 


r din Zg 

-9 V 


dfx 


( 18 ) 


B. The Traceless Part of the Energy-Momentum Tensor 

The traceless part of the energy momentum tensor can be written as a sum of four parts, 


rnjilf _ rp[ii/ , rrqus , rp/iv , 

1 — 1 g \ 1 q ' 1 gv ' 


( 22 ) 


The gauge-invariant gluon part Tjf u is, 


j? u ) 

± a L L a i 


(23) 


where and henceforth the symbol (////) also means subtraction of the trace. The gauge- 
invariant quark part is 




( 24 ) 
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The gauge variant part is, 


Tg = C l ^d u) d • A + d (fi idD v ^uj , (25) 

which is the BRST variation of d^CuA u \ Finally, the gluon-equation-of-motion-related op¬ 
erator is, 


E^ 


- A 


6S 


(26) 


which is finite. 

According to 
renormalization. 


the Joglekar and Lee theorems [16], the above four operators close under 


Furthermore, the mixing matrix takes the following form, 


( T r\ 


( Zqq 

Zqg 

Zqa 

Zqe ^ 


( T r\ 

g 


Zgq 

Zgg 

Zga 

Zge 


T^ v 

g 

gv 


0 

0 

Z a a 

Z a e 


gV 

\E^ ) 


l o 

0 

0 

1 ) 


V E> lv / 


(27) 


That is, the gauge variant operators do not mix with gauge invariant operators and the 
equations of motion operators do not mix with any other operators. 

To find the Z’s, one needs to study Green’s functions with one insertion of the composite 
operators in perturbation theory. However, we can find constraints among the renormal¬ 
ization constants by studying Ward identities related to space-time translational symmetry. 
Calculating the divergence of the energy-momentum tensor without using the equations of 
motion, we find, 


= 


SS a „ AU SS 5S auJ SS 5S^ U 

d A >1 - —d v if; - -jd''ip - —d v u - —d v u) . 


5A 


dtp 


Sip 


Scj 


5cu 


(28) 


When the above equation is inserted into a renormalized Green’s function (Ward identity), 
the right-hand side simply replaces the elementary fields by their derivatives, and thus is 
finite. Hence the divergence of the energy momentum tensor is a finite operator. The only 
remaining counter term for the tensor itself is of form 

(S"S“ -sr&)f {A,U,U>) , (29) 


where / is a dimension-two function of its arguments. However, since the trace part of 
the energy-momentum tensor is finite, Lorentz symmetry allows only counter-term tensors 
that are symmetric and traceless. Thus the traceless part of the tensor must be finite by 

itself 0. 

The finite traceless part of the energy-momentum tensor imposes the following con¬ 
straints on the renormalization constants in minimal subtraction scheme, 


Zqq T Zgq f ■ 

Z qg + Z gg = 1 , 
Z qa + Z ga + Z aa = 1 , 
Z qe + Z ge + Z ae = 0 . 


(30) 
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Thus in the scheme, we have, 

T£ v + + Tg + E^ = T^ R (g 2 ) + T^ r (p 2 ) + T^ r (/i 2 ) + E^ . (31) 

Although individual term on the right hand side depends on the renormalization scale, the 
sum does not. Again, according the Joglekar-Lee theorems, Tj^ R (g 2 ) and E^ u have vanishing 
physical matrix elements, and for practical purposes, we can keep only the gauge-invariant 
quark and gluon contributions in the traceless part of the energy-momentum tensor, 


T'"' = T'“' R (fi 2 ) + T^ R {p?) . 


(32) 


Recently, some questions arise in the literature about validity of the Joglekar and Lee the¬ 
orems [[T0l.[20|| , in particular, regarding the form of the mixing matrix appearing in Eq. (^7|) . 
Harris and Smith [pTI] pointed out that so long as one is working with composite operators at 


non-zero momentum transfer, all the theorems should remain valid. Collins and Scalise [20 


on the other hand, have worked at zero momentum transfor. They showed that the on-shell 
limit for gauge bosons is subtle and potentially causes problems. However, it appears that 
if one works with off-shell Green’s functions or physical hadron states, the operator mixing 
shall follow the standard theorems. More work is certainly needed in this direction to clear 
up the issue. 


C. Scale Symmetry and the Anomalous Ward Identities 


The above discussion on mixing of dimension-4 operators and the trace of the energy mo¬ 
mentum tensor permits a simple explanation for the anomalous breaking of scale symmetry 
and a derivation of Callan-Symanzik equation (or the anomalous Ward identity). Consider 
the scale (dilatation) transformation in d dimension, 


x - 

-> Ax , 


i>{x) - 

d-1 

-> A 2 7p(Ax) , 


A^(x) - 

-4 A^A^Ax) , 


u;(x) - 

A^ -1 u;(Aa;) . 

(33) 


The change in the effective QCD lagrangian is a total derivative plus symmetry breaking 
terms. There are two types of symmetry breaking effects: quark masses and the gauge 
coupling, the latter has dimension e/2 in dimensional regularization. Thus, 


e BE 

5E eS = d ll (x tl E eS ) + - -g—^- • (34) 

2 dg 

The last term is clearly not gauge invariant. A simple rearrangement reveals that it consists 
of operators Op and O g f defined in the previous section plus a total derivative term, 


~9 


dCcft 

dg 


d 11 (F^A 1 ' + C\d • A)A„) + 2 0 F + 2O gf . 


(35) 
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Since 0 F and O g f diverges like 1/e, the scale symmetry is broken anomalously at the quantum 
level. 

On the other hand, with use of the equations of motion, one has, 


= (36) 

where (f> t is a generic notation for all the fields. According to this, we can define a dilatation 
current corresponding to the scale transformation, 

J d = ^ (F^A„ + • A) A*) . (37) 

odpQi 2 v J 

The symmetry breaking terms now appear as divergence of the current, 

d^F, = %pimp + e (O f + Ogf) . (38) 

To prove that, use of the equations of motion is essential. 

The dilatation current has a simple relation with the energy-momentum tensor in Bclin- 
fante’s form. To see that, we use the definition of the canonical energy-momentum tensor 
Tq' to write, 

Jg = TFx u + gAj-(>"/!„ + r'(3 ■ A)A») , (39) 

where d$ is the canonical dimension for pi field. Using the Belinfante improvement, one has, 

r n = t^x„ + + <r<4),A - - (f^a* + r\a■ a)a„) , («) 

where the second term is called held-virial, which has a part that cancels the last gauge 
non-invariant term, 

= T^x u + (i - |) (d^uxjj + ujD^uj) + (2 - e)C\d ■ A)A fl . (41) 

The remaining terms cannot be written as a total derivative. According to Callan, Coleman, 
and Jackiw 0, this means that conformal symmetry present in the basic QCD lagrangian is 
broken by the gauge fixing. This fact is of course well-known in perturbation theory where 
a gauge-fixed gluon propagator does not have conformal symmetry. However, the breaking 
term is a BRST-exact operator plus the ghost equation of motion and thus has no physical 
consequences. 

To derive the anomalous Ward identities associated with scale transformation, we con¬ 
sider corresponding change in the Green’s function generating functional Z( J M , 77 , 77 , x, X)- A 
simple derivation yields the following equation, 



G R (xj) + G r ( Xj ,&J Dli ) 


0, 


(42) 
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where G R (xi ) is a renormalized Green’s function with n a external gluon lines, quark 
(antiquark) lines, and n w ghost (antighost) lines. G R (xi, is the same Green’s function 

inserted with divergence of the dilatation current. In momentum space, we have, 


- 4 ( n ^ + n^ - 1) + (n A + | n^ + n u ) 

■ opi z 


G R (pi) + G u (pi, — 0 . (43) 


-iR 


Here momenta pi are conjugation of all Xi except for one which is taken to be zero. On the 
other hand, a simple dimensional analysis yields, 


d d d 3 

E Pi-fy. + ^ + mR Q^R + + n P + ( rlA + 2 n * + n ") 


G R (Pi ) = 0 , 


(44) 


where p is a renormalization scale. Combining the above two equations, we get, 


d 


,R 


d 


% + m dm R 


G R (p i ) + +G K {p„d>‘J Dr ) = 0. 


iR 


(45) 


According to Ref. |TT[ , an insertion of the divergence of the dilatation current can be re¬ 
placed by derivatives with respect to the gauge coupling, quark masses, and the gauge-fixing 
parameter. Thus, we find the following anomalous ward identity (Callan-Symanzik equa¬ 
tion), 


d d 

^ +>} dpt 


- 7 m'm 


R 


d 


dm R 


-7s£ 


R 


d , 7s . 72 . 7 

dt* +nA Y + n *Y + n “2 


G R (pi) = 0 . (46) 


This is just the well-known renormalization group equation which can be derived indepen¬ 
dently by studying the dependence of the renormalized Green’s function on the renormal¬ 
ization scale. 


III. BREAKUP OF THE HADRON MASSES 

Let us first recapitulate the main results obtained in the last section. First of all, we have 
shown explicitly that the QCD energy momentum tensor T^ u (with vacuum subtraction) 
is a finite composite operator and thus a physical observable. Second, the tensor can be 
separated uniquely into the traceless and trace parts, each of which is a finite and scale- 
independent composite operator, 


rjifll' _ JiflU _|_ Ji/ih 


(47) 


So the separation at this level is completely physical. Finally, in physical states the traceless 
part of the tensor is effectively a sum of the renormalized, gauge-invariant quark and gluon 
contributions, 


= T^ R (p 2 ) + T^ R (p 2 ) , 


(48) 
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and the trace part is a sum of the gauge-invariant quark mass and trace anomaly contribu¬ 
tions, 

= T^ R (ii 2 ) + T^(/i 2 ) . (49) 

The separation at this second level is renormalization and regularization scheme dependent. 

In this section, we study a breakup of hadron masses according to the above properties of 
the energy-momentum tensor. We first discuss matrix elements of various parts of the tensor 
in hadron states. Then we write down the finite form of the QCD hamiltonian and drive a 
mass separation formula. We revisit its application to the nucleon and discuss implications 
of the result on the nucleon’s quark-gluon substructure and on the low-energy dynamics of 
QCD. Application to the pion is studied in a separated subsection, where we show that the 
color electric and magnetic fields in the pion is the same as those in the QCD vacuum in 
chiral limit. We also discuss the color fields in the QCD vacuum in a Lorentz covariant, 
non-perturbative regularization scheme. 


A. Matrix Element of the Energy-Momentum Tensor in Hadrons 

Let us consider the forward matrix element of the energy-momentum tensor in a hadron 
state, \P ), where P^ is the four-momentum of the state. We assume the state is normalized 
according to ( P\P) = (E/M) (27t) 3 <5 3 (0), where E = P° is the energy of the state and M is 
the mass of the hadron. It is well known that, 

(P\T^\P) = P^P U /M . (50) 

A simple derivation of the above equation goes like this: According to Lorentz symmetry, 
one has, 


(P\T^\P) = aP fl P l/ + bg^ . 


(51) 


where a and b are scalar constants. On the other hand, the hamiltonian of the system is 
H = f d 3 x T 00 , which has the following matrix element in the hadron state, 

(P\H\P) = (E 2 /M)(2n) 3 5 3 (0) . (52) 


Comparing Eq. ( |52| ) with Eq. we obtain Eq. (|50|) . 

Eq. m) allows us to obtain the matrix elements of the traceless and trace parts of the 
energy-momentum tensor separately. In fact, decomposing both sides of Eq. (j5^) into trace 
and traceless parts, we have, 


(P\T^\P) = (P»P U - /M , 

(P\T^ U \P) = . (53) 


The right-hand sides of both equations are independent of the renormalization scale, as 
required by Lorentz symmetry. The second line in (|53|) is not new; it has appeared before 


in the literature 22-24 
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We use Lorentz symmetry again to define two scheme-dependent matrix elements below. 
First, we define matrix element of the quark operator appearing in the traceless tensor, 

(P\T^\P) = a(n 2 ) (P^P U - ^M 2 ) JM , (54) 

where a(/i 2 ) has an explicit scale-dependence. From Eq. (|48D and the hrst line in Eq. (p3|), 
we have, 


(/'|7J‘" /') = (l - ain 1 )) - ijC-M 2 ) /M . 


(55) 


Second, we dehne matrix element of the quark-mass operator, 


(p|frip> = b{^)-sTM , 


(56) 


where the renormalization scale-dependence comes entirely from the anomalous dimension 
7 m of the mass operator, which depends on the renormalized gauge coupling. From Eq. 
and the second line in (|53|), we have, 


<p\fr\p) = (i - f>(/7) Vm. 


(57) 


Thus the matrix elements of the four parts of the energy-momentum tensor are entirely 
determined by two parameters a(/P) and 5(/i 2 ). 

The matrix element a(/i 2 ) can be extracted from hadron structure functions measured in 
deep inelastic scattering. According to operator product expansion, the twist-two operator 
Tq U appears in the expansion of product of two vector or axial-vector currents. Using a 
dispersion relation, one can relate the matrix element of T£ v to the hrst moment of quark 
distributions, 


a(/i 2 


= £ 

/ ■ 


dxx 


q/(x, 7) + qf(x, y 


(58) 


where qf and qf are quark and antiquark distributions in the hadron. The physical meaning 
of a(/i 2 ) is the momentum fraction of the hadron carried by quarks in the infinite momentum 
frame |25j . 

The matrix element b(/i 2 ) is proportional to the nucleon’s scalar charge (P|^^|P), which 
cannot be measured directly from an experiment. However, the operator is a part of the 
QCD lagrangian which breaks the chiral symmetry explicitly. Therefore, useful information 
about b(/j, 2 ) may be obtained by studying physical effects of chiral symmetry breaking. 


B. A Separation of Hadron Masses 
The QCD hamiltonian is defined as, 

Pqcd = J d 3 xT 00 (0,x). (59) 
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From the discussion of the last section, Hqco is a finite operator. [Interestingly, however, 
the lagrangian density has 1/e divergences.] According to the separation of the energy- 
momentum tensor, we have the partition of the hamiltonian operator, 


Hqcd — H' q + H g + H' m + H a , 

where various terms are, 


K 

H g 

H a 



3 - 

ijv{—iD ■ a)^> + -ilmu/j 



j d 3 x-(\ + r y m )ipmip , 
I d 3 x i/3( S )(E 2 -B 2 ). 


Rearranging the mass term by defining, 


H q = J d?x^\—iD ■ a)ij} , 

H m = j d?x (l + ipmip , 


(60) 


(61) 


(62) 


we have Eq. (|60[) without the primes. Note that the above expression implicitly contains 
the renormalization counter terms, which we omit for simplicity. The physical meaning of 
the various pieces of the hamiltonian is clear: H q represents the quark kinetic plus potential 
energy (the static color interaction is not included); H g is the gluon kinetic and potential 
energy; H m is the quark mass contribution; and finally, H a is the anomaly contribution 
whose significance will become clear later. 

In field theory, masses of bound states are usually defined as poles in Green’s func¬ 
tions (implicitly so in lattice calculations). In hamiltonian formalism, one can define masses 
as eigenvalues of a hamiltonian as in non-relativistic quantum mechanics. In the follow¬ 
ing discussion, we assume that hadron masses are calculated as expectation values of the 
hamiltonian operator in the rest frame of hadrons, 


M 


(m Q cp|p) 

(p\p) 


rest frame • 


(63) 


We view the above equation as a mass probe into the structure of quark and gluon config¬ 
urations in hadron states, since different pieces of the hamiltonian have different sensitivity 
to various components of hadron wave functions. At this point, the reader might recall the 
recent investigation in the literature about the spin structure of the nucleon |22|] , where one 
is interested in how the spin of the nucleon is made of the spin and orbital angular momenta 
of quarks and gluons. The same question can be asked of the mass structure of hadrons. 

According to Eqs. (0), (0, and (163[) , masses of hadrons can be separated into various 
contributions, 


M — Mq + Adg + M m + M a . 


(64) 
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Relating the matrix elements of different pieces of the hamiltonian to those of the energy- 
momentum tensor, we have, 


M q 

M g 

M m 

M a 


4 \ 1 + 7 m ) 


4 + 7m 

4(1 + 7m) 


bM , 


-(1 — b)M . 


M, 


(65) 


Therefore, knowing the parameters a and 6, we can determine a complete separation of the 
masses. In the following subsections, we shall apply the separation to several cases and 
discuss physical significance of the result. 


C. The Mass Structure of the Nucleon 


The mass structure of the nucleon was analyzed in detail in the Letter paper 0. Let us 
briefly summarize the main result here and discuss its physical implications more thoroughly. 
In the end, we shall make some remarks for the mass separation of other hadrons following 
a rather plausible assumption. 

The matrix element a(/i 2 ) for the nucleon has been measured quite accurately in various 
high-energy scattering processes involving nucleons. A recent extraction gives 


*MS 


(1 GeV 2 ) = 0.55 . 


( 66 ) 


In the Letter, two estimates of 5(1 GeV 2 ), in the limits of chiral SU(3) ( m s —> 0) and heavy 
strange quark (m s —> oo), were given. They are, 


b(m s -> 0) = 0.17 , 

b(m s —> oo) = 0.11 . (67) 


Since the two limits do not lead to qualitatively different conclusions, we focus on the chiral 
limit below. The numerical numbers in this case are, 


M q = 270 MeV , 

M m = 160 MeV , 

M g = 320 MeV , 

M a = 190 MeV . (68) 

According to the above result, the quark kinetic and potential energies contribute about 
1/3 of the nucleon mass. The further separation into the two is not gauge invariant. However, 
the practice does have an intuitive appealing at the phenomenological level. If there are three 
massless quarks confined to a spherical cavity of radius 1 fm, as in the MIT bag model, the 
total kinetic energy is 600 ~ 700 MeV. Thus, the color-current interaction between quarks 
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and gauge fields contributes —300 ~ —400 MeV to the mass, which is consistent with the 
magnitude of N — A splitting. Such large current interaction is intrinsic to a truly relativistic 
theory. It does not occur for instance in low-energy QED, where the static Coulomb potential 
plays a dominant role. The strong current interaction certainly induces quark interactions 
of Nambu-Jona-Lasino type, and is perhaps at the origin of the chiral symmetry breaking. 

The quark energy can be further separated into contributions of different flavors. The 
parameter a measures the fraction of the nucleon momentum carried by quarks, which is 
known separately for each flavor. For instance, 0.38 fraction of the nucleon momentum is 
carried by up quarks, which can be translated into the up quark contribution to the nucleon 
mass 250 MeV. Likewise, down quarks contribute 105 MeV mass, and strange quark pairs 
contribute —85 MeV. One might try to break the contribution from each flavor into these 
from valence and sea quarks. However, since the valence and sea contributions to the matrix 
element b are unknown, the separation is not possible. Nonetheless, for up and down quarks, 
the contribution from b is not large, and may be neglected. From the momentum fractions 
carried by the sea, we find that the up or down sea contribution to the nucleon mass is on the 
level of 30 MeV. The small number indicates that there are not many quark and antiquark 
pairs in the wave function. Thus the nucleon seems to have a simpler Fock expansion than 
the small current quark masses imply. 

The quark mass term contributes about 1/8 of the nucleon mass. About half of which 
comes from the strange quark pairs. The strange quark contribution here is definitely less 
certain than two light flavors. One might hope that future lattice measurement of the strange 
scalar charge may reduce the uncertainty. The implied magnitude of the scalar charges 
(P\uu\P), (P\dd\P), which are charge conjugation odd quantities, is another indication that 
the number of quark-antiquark pairs is small. 

The normal gluon energy contributes about 1/3 of the nucleon mass. This energy includes 
the color-static Coulomb energy between quarks. The gluon part of the trace anomaly 
contributes about 1/4 of the mass. From these, we deduce the color-electric and color- 
magnetic fields in the nucleon separately (taking a s (l GeV 2 ) = 0.4), 

(P|E 2 |P) = 1700 MeV , 

(P|B 2 |P) = -1050 MeV . (69) 

The second line indicates that the color magnetic held in the nucleon is smaller than that in 
the vacuum. This property of the magnetic held has long been suspected phenomenologically. 
The present result lends a strong support for the educated guess. Clearly, this behavior of 
color helds is closely related to color confinement. 

In the chiral limit, the trace anomaly contribution is analogous to the vacuum energy in 
the MIT bag model. In fact, the trace part of the energy-momentum tensor in the bag is 
Bg MI/ , where B is the energy density of the “perturbative vacuum.” The role of such energy 
density is to confine quarks. Thus the scale symmetry breaking is explicitly connected 
to quark confinement. It is essential then to include the effects of the trace anomaly in 
phenomenological hadron models. 

A final comment is about the role of strange quarks. They contribute —60 MeV to the 
mass through the trace anomaly because the [3 function depends on the number of flavors. 
The kinetic and potential energy contributes about —85 MeV. Adding these to the strange 
mass contribution of 115 MeV, one gets a total of —30 MeV, roughly three percent of the 
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total mass. Therefore the uncertainty in the separation is largely limited to the strange 
sector. 

What shall be the general feature of the mass separation for other hadrons for which 
there are no data? First of all, so long as one is concerned with non-strange hadrons, the 
contribution of the quark mass term is presumably small and we may neglect b. Second, that 
gluons carry about half of the nucleon momentum in infinite momentum frame is perhaps 
approximately true for all hadrons, e.g., p or A. Thus we further assume a = 0.5. Given 
these guesses, we have, 

M q = ^M, M g = hd , M a = X -M . (70) 

This is a heuristic way to sum up the main result of the mass separation. 


D. The Mass Structure of the Pion and the Vacuum 


The mass structure of the pion is particularly interesting because, according to the 
Goldstone theorem, the pion is intrinsically different from ordinary hadrons: it is a collective 
mode in the QCD vacuum. As we shall see, the mass structure indeed reflects this. 

The matrix element a(/i 2 ) can be extracted from the quark distributions measured in 
7 t — N Drell-Yan processes. The quality of the available data |27|, however, is much less 
satisfactory compared with that of the quark distributions in the nucleon. Nonetheless, it 
seems safe to conclude that a(l GeV 2 ) is known at ten percent level, with a central value 
similar to that of the nucleon, 


a(l GeV 2 ) = 0.55 ± 0.05 . (71) 

[Note that the precision of the data is not good enough to discern radiative corrections at 
the sub-leading-logarithmic level, so we have suppressed the scheme label.] Thus as far as 
high energy probes are concerned, the pion is not dramatically different from other hadrons. 
This is also true for the i t — N total cross section at high energy, for which the quark counting 
rule appears valid. 

The matrix element b can be calculated through study of the pion mass in chiral pertur¬ 
bation theory. To avoid kinematic singularity in chiral limit, we adopt the normalization of 
the 7r state, (P\P) = 2P(27t) 3 <5 3 (0). Therefore the matrix element b is, 


bm 2 n = (P | ^ m f ^ f ^ f \P 

\ f=u,d,s / 


On other other hand, the first-order chiral perturbation theory predicts |25 

ml = ~{m u + m d )(0\uu + dd\0)/fl . 


(72) 


(73) 


where /k is the pion decay constant and |0) is the QCD vacuum. The above equation tells 
us two things. First strange quarks do not contribute to the pion mass in the first-order 
perturbation, and matrix element (P|ss|P) is strongly suppressed. Second, m 2 ~ m u ,md, 
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and so to the first order accuracy, one can take the chiral limit of the pion wave function in 
the right hand side of Eq. (72). 

The pion mass can also be calculated by using the ordinary first-order perturbation 
theory p8| , 


-ml- = [P\m u uu + mddd\P^j , 


(74) 


where | P) is the pion wave function in chiral limit. Comparing the above equation with 
Eq. (|72|), we have, 


6 = 1 / 2 . 


(75) 


So the first-order chiral perturbation gives a clean prediction. 
Using the above matrix elements, we roughly have, 


M q — 0 , 


M g = 




47 m 2 m ^ 


Mn = 




(76) 


Two comments can be made immediately with regard to the above mass partition. First, the 
quark kinetic and potential energies cancel almost exactly. This fact is difficult to reproduce 
in quark models for the pion, where quarks carry large kinetic energy when confined to a 
small region of space. Second, the color electric and magnetic fields in the pion approach 
those in the vacuum in chiral limit. This strongly indicates that the pion is a collective 
excitation of the QCD vacuum. 

It is tempting to use the above formalism to study the color electric and magnetic fields 
in the QCD vacuum. However, the energy-momentum tensor without vacuum subtraction 
(“normal-ordering”) is not finite and it is dangerous to work with divergent quantities. 
Nonetheless, if one finds a Lorentz covariant and non-perturbative regularization scheme to 
define the energy-momentum tensor without the subtraction, one can make statements about 
the vacuum color fields, except, of course, discussion is regularization scheme-dependent and 
is meaningful only in a carefully defined context. [Perturbative regularization schemes are 
usually plagued with the infrared renormalon problem, and hence are not useful in this 
context 0] 

For simplicity, let us neglect quarks and study pure non-Abelian gauge theory. Since the 
vacuum is not characterized by any four-vector, the matrix element of the traceless part of 
the energy momentum tensor vanishes. That is, 


(0|T^|0) = 0 . 


(77) 


Taking /j, = u = 0, one gets, 


(01E 2 10) = — (01B 2 10) . (78) 

The color electric held is negative of the color magnetic held in vacuum! This appears to 
be a quite dramatic statement at hrst. However, it is trivially true, with some caveats, in 
lattice QCD calculation . On lattice, the color electric and magnetic helds can be defined as 
the average of the trace of elementary plaquettes in space-time and space-space planes. Due 
to hypercubic symmetry on lattice, the two different plaquettes have the same expectation 
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value in vacuum. Remembering that the electric field in Minkowski space is related to 
that in Euclidean space by a factor of i, the imaginary unit, one gets the above equation 
immediately. 

The trace part of the energy-moment urn tensor contains only the anomaly. According 
to Lorentz symmetry, 


(0|f^|0) = g^p 


(79) 


where p is the vacuum energy-density in a particular regularization scheme. Using the trace 
anomaly, one has, 

0|(-2/3/^)F^F^|0) = 4p (80) 

The above equation shows that the vacuum energy density is related to the expectation 
value of F 2 , which is called the gluon condensate in the literature. Knowing the condensate 
in a particular scheme, one obtains the color electric and magnetic fields of the vacuum, and 
hence the vacuum energy density. 

Unfortunately, no one has yet proposed a natural Lorentz covariant, non-perturbative 
regularization scheme for the unsubtracted energy-momentum tensor operator, ffowever, 
there is a phenomenological definition of the vacuum condensate used in the QCD sum rule 
calculation H. The magnitude of the condensate has been determined through fitting to 
hadron spectra, 


(o|^F 2 |o\ = (0.35 GeV) 4 . 


(81) 


This translates to, 


(0|B 2 |0) = —<0|E 2 |0) = 3.68 GeV/frn 3 . (82) 

It is a rather large number from the phenomenological point of view. However, one should 
not forget its definition when it is used in a physical context. 


IV. DISCUSSIONS AND COMMENTS 

The structure of the nucleon is a subject that has been discussed for many years. The 
non-relativistic quark model has the virtue that it is simple and captures many important 
aspects of physics. Unfortunately, to improve our understanding, we must solve QCD in 
non-perturbative region. Although lattice QCD provides an effective method to calculate 
many observables, it provides little insight about the physics. 

The measurement of the spin structure functions of the nucleon is a milestone in motivat¬ 
ing new explorations of the quark-gluon structure of the nucleon. It points to our deficiency 
in traditional modeling of hadrons and to a need for unquenched lattice calculations. On 
the other hand, it also urges a better physical description for hadrons. In light of this, any 
rigorous information about the nucleon properties is useful. 

The mass separation of the nucleon is certainly one step towards a better understanding 
the quark and gluon dynamics of the nucleon. As was discussed, it has many implications 
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about the physics of gluons and quarks, and their interactions. What is particularly inter¬ 
esting is the anomaly contribution. If the reader is familar with the “spin crisis,” he/she 
might recall that the axial anomaly was considered as one of the constributions to the nu¬ 
cleon’s spin |H|. Unfortunately, there the separation between the anomaly and normal 


contributions are not quite clear. In particular, the question of gauge invariance and factor¬ 
ization has not be solved satisfactorily. Here the anomaly contribution to the nucleon mass 
is unambiguous and its physics interpretation is quite clear. 

One may ask where one goes from here. First, one can confirm the present result by 
doing lattice calculations. One can measure different pieces of the hamiltonian in the nucleon 
state. Such a calculation in the end may help to test lattice approximations. Second, one can 
try to build models which are consistent with the present mass separation. One important 
conclusion here is that the anomaly term must be added to models like Nambu-Jona-Lasino. 
Without this term, it is difficult to include the confinement effects. 
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